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Abstract. 

We theoretically analyze equilibrium fluctuations of persistent current (PC) in 
nanorings. We demonstrate that these fluctuations persist down to zero temperature 
provided the current operator does not commute with the total Hamiltonian of the 
system. For a model of a quantum particle on a ring we explicitly evaluate PC noise 
power which has the form of sharp peaks at frequencies set by the corresponding 
interlevel distances. In rings with many conducting channels a much smoother and 
broader PC noise spectrum is expected. A specific feature of PC noise is that its 
spectrum can be tuned by an external magnetic flux indicating the presence of quantum 
coherence in the system. 
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1. Introduction 

Meso- and nanorings formed by normal conductors and pierced by external magnetic 
flux develop persistent currents |1]. This fundamentally important equilibrium effect 
is a direct consequence of quantum coherence of electrons which - at low enough 
temperatures - may persist up to distances exceeding the perimeter of such rings. 

Does persistent current (PC) fluctuate? At the first sight it might appear reasonable 
to assume that at least for T — j- no such fluctuations could occur. Indeed, while at non- 
zero T thermal fluctuations of PC should be present [2], in the zero temperature limit 
the system approaches its (non-degenerate) ground state and, hence, no PC fluctuations 
would be possible. 

Below we will demonstrate that in many cases it is not so. Namely, no PC 
fluctuations are expected in the zero temperature limit only provided the current 
operator commutes with the total Hamiltonian of the system, otherwise fluctuations 
of persistent current can occur even in the ground state exactly at T = 0. Theoretical 
and experimental investigations of such PC fluctuations can give important additional 
information about the ground state properties of meso- and nanorings not contained in 
the average value of PC. 

Note that fluctuations of PC in the ground state may be induced provided the ring 
interacts with some quantum dissipative environment. In this case such interaction 
is responsible for (i) phase-breaking effects implying suppression of both quantum 
coherence and PC and (ii) non-vanishing fluctuations of PC down to T — )■ 0. E.g. 
it was demonstrated p] that interaction with Caldeira-Leggett environment decreases 
the average value of PC and simultaneously increases PC fluctuations which are directly 
related to fluctuations in the environment itself. Such effects are of importance, e.g., 
for PC qubits which quantum states can be entangled with those of environment [1]. 
Fluctuations of PC down to T — > could also occur provided the number of particles 
in a ring fluctuates due to its interaction with some reservoir jEj- 

The situation considered here is entirely different: We do not assume the presence 
of interaction or particle exchange with any environment at all. Accordingly, quantum 
coherence of the system is fully preserved and no PC suppression takes place. As it 
will be demonstrated below, quantum coherence implies the possibility of tuning of PC 
fluctuations by an external magnetic flux applied to the system. 

2. The model and general relations 

For definiteness, let us consider a simple model of a quantum particle with mass M on a 
Id ring of radius R pierced by magnetic flux $, see, e.g.. Fig. 1. The particle position on 
the ring is parametrized by the angle 6 which will be the quantum mechanical variable 
of interest in our problem. The Hamiltonian for this system reads 




(1) 
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Figure 1. The system under consideration: A particle on a ring in the presence of a 
periodic potential. The ring is pierced by the magnetic flux. 



where = — is the (dimensionless) flux operator, U{6) defines the potential profile 
for our particle, (px = ^/^o and $o is the fiux quantum. Within this model PC was 
previously studied in a number of papers P, [3, El El [ini [H] in the presence of various 
dissipative environments. In addition, the model discussed here could be of interest for 
the problem of PC in superconducting nanorings in the presence of quantum phase slips 

m- 

For our problem the current operator in the standard Schrodinger representation 
is defined as 

Switching to the Matsubara representation 

/M(r) = e-^/e-^^, (3) 
we define the current-current correlation function 

oo 

n(r) = ((r/M(r)/M(o))) = T y: n,^„e— , (4) 

?i=— oo 

which describes equilibrium current noise. Here T is the time-ordering operator, 
Un = 27mT are Matsubara frequencies, (...) = tr(p...) denotes averaging with the 
equilibrium density matrix p = e~^^/Z, where Z = tie~^^ is the grand partition 
function and /3 = 1/T. The symbol ((. . .)) stands for irreducible correlators (cumulants), 

e.g., ((/(r)/(o))) = mm) - {f{T)){m) . 

Employing the full set of eigenstates H\m) = after a straightforward 

calculation we obtain 

n(r)=P + n(r), (5) 

where V does not depend on imaginary time and reads 

p = l^|(m|/|m)|V^^'" - (6) 



Fluctuations of PC 



4 



while the Fourier components of n(r) are defined as 

In order to estabhsh the relation between the correlator (jl]) and the current noise 
power we also define the Heisenberg operators /(t) = e**^Je~**''^ and the Keldysh Green 
function 

S{t) = ((/(t)/(0) + /(0)/(t))) = / ^5.e-*. (8) 

As before, decomposing the result for S{t) into time-independent and time-dependent 
contributions we find 

Sit) = 2V + Sit), (9) 

where V is again defined in Eq. ([6]) and the Fourier components of S{t) take the form 

^. = ^E iH/in)r 

X (e-^^- + e-^'-) 6{io + - (10) 
Comparing now Eqs. (j7]) and (ITUI) we arrive at the fiuctuation-dissipation relation 

S^ = 2 coth ^^^n^^+io, (11) 

which allows to immediately recover the current noise power from the correlator 
Note that this relation links together the quantities 11 and 5* which - according to Eqs. 
(jS]) and dH]) - differ from the the correlators 11 and S by the constant in time terms, 
respectively V and 2V, which produce singularities in the frequency domain. 

We would like to point out that our formalism also allows to analyze the linear 
current response to the time- dependent fiux inside the ring and to formally define the 
ac conductance of the system. According to the Kubo formula this ac conductance 
is expressed in terms of the commutator of the current operators, unlike the noise 
spectrum defined by the anticommutator of these operators ([8]). Below we restrict our 
attention only to time-independent values (p^ and do not address the behavior of the ac 
conductance which we are not interested in here. 

The above exact relations fully determine PC correlators in terms of the system 
eigenstates. These relations allow to observe that as long as the current operator / 
commutes with the Hamiltonian of the system the Fourier components ([7]) and ffTOl) 
vanish identically together with the matrix elements {m\I\n) with m ^ n, while the 
time- independent term V (j6]) tends to zero only in the zero temperature limit. Thus, 
in this case no PC fiuctuations can occur at T — )■ and at non-zero temperatures PC 
noise does not vanish only in the zero frequency limit, = 2V6{uj). 

If, however, the operators I and H do not commute with each other the situation 
becomes entirely different. In that case the matrix elements {m\I\n) in general remain 
non-zero for any pair of eigenstates and, hence, PC fluctuations may persist down to 
T = 0. 
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For a simple model of Fig. 1 PC correlators can be evaluated directly from Eqs. (jl])- 
ffTOj) . In more complicated situations, however, the above general relations employing 
the representation of eigenstates could become less convenient for practical calculations. 
For this reason below we will develop alternative approaches which can also be useful 
for the analysis of PC fluctuations. 



3. Free energy and current noise 

It turns out that in both limits of zero Matsubara frequency and zero imaginary time the 
current-current correlator can be conveniently related to the free energy of the system 
J-' = —TlnZ. Making use of the expression for PC 

(/) = -1^ (12) 
together with the identity 

^ = _/,„-'-"*|^e-* (13) 



for the second derivative of the free energy with respect to the flux we obtain 

- ^ ~ dr^ ^ ^ + /?(/)'. (14) 



Att^ d(f)l AttMR^ J tiie-'^H 
From this equation one readily finds 



Ati^MR^ 47r2 d(t)l 

u 

On the other hand, employing the identity 



we get 



and, hence, 



= (17) 

As in many cases the free energy of the system can be readily evaluated, Eqs. f|T5|) 
and flTSl) provide a great deal of information about PC noise. Further simplifications 
may appear in the zero temperature limit since in this case the free energy reduces to 
the ground state energy F{T — )■ 0) = £o(0x)- E-g- for a free particle on a ring (i.e. for 
U (6) = 0) one has £o(</'x) = 0^/ {2MR'^) and, hence, in this case in the limit T ^ from 
Eqs. (|T5|) and f lTSj) one trivially finds 

Ho = n(0) = 0. (19) 
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In agreement with our general analysis in the absence of an external potential and 
at T = the correlator (jl]) vanishes identically for all values of r implying that no 
fluctuations of PC occur in this case. This is because for U{9) = the current operator 
commutes with the Hamiltonian. 

At non-zero external potentials U{6) ^ 0, however, these two operators do not 
commute anymore and, hence, fluctuations of PC in general do not vanish even at very 
low T. This conclusion can be reached, e.g., from Eq. f|T5|) without any additional 
calculation. Indeed, for U{9) ^ the ground state energy eo(0x) deviates from 
(l)^/{2MR^) and, hence, Hq ^ down to T = 0. 



4. Generating functional 

Let us now formulate a general technique that will allow to fully describe current 
fluctuations of PC within our model. For this purpose we define the generating 
functional 

Z[7]]= V(t)Veei ^ \ (20) 



where r]{T) is the source field for the flux variable (p. Performing integration over (p we 
obtain 

Z[r]] ^ Vee ^ ' . (21) 



Taking the variational derivative of J-'[?7] = —T\YiZ\!]\ over the source field ?7(r) and 
setting this field equal to zero afterwards, we derive the relation between the expectation 
values for the current and the particle "velocity" Q: 

IP 

(Hr)) = ^(^(r)). (22) 

Similarly, the second derivative of J-'[r]] with respect to ?7(r) yields the second current 
cumulant: 

-^mri)0(T,m. (23) 

Analogously one can establish the relations between higher current and velocity 
cumulants. Up to some unimportant (5-functions at coinciding times (which cancel out 
in the final result as we will see below) the latter cumulants, in turn, are evaluated from 
the relation 

mM.m = (-if^£^i§^^u^_„, (24) 

where Z[Q is the generating functional 

Z[Q= / dOo J2 e'"""'^^ (25) 
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The above general expressions allow for straightforward evaluation of all current 
cumulants thus establishing "full- counting statistics" of PC in our problem. 

5. Current-current correlator 

Below we will focus our attention on the current- current correlation function (jl]) which 
will be evaluated in the specific limiting case [11] 

U{e) = Uo{l - cos(Ke)), Uo > K^/{MR^). (26) 

In other words, we will assume that the particle confined to the ID ring is moving in 
a periodic potential with the distance 2ti / n between adjacent minima. For k = 1 our 
model reduces to that derived for ultra-thin superconducting rings in the presence of 
quantum phase slips |12]. As indicated in Eq. fl26l) the potential barriers between these 
minima are high, in which case the particle moves around the ring due to hopping from 
one minimum to another. Semiclassically, these hops are described by multi-instanton 
trajectories [11] 

Q{r) = Y.vfe{r-T,), z/, = ±l, (27) 

i 

which dominate the path integral (!26]l . Here 9{t) = 4 arctan(e^'^)/K is well known kink 
solution, describing the particle tunneling with the amplitude 

A/2 = 4(r]f/o/7r)i/2e-^, (28) 



where f2 = k^JUq/{MB?). Substituting the trajectories ( 1271) into (l26l) and performing 
Gaussian integration we get 

/3 /3 



n=0 ui=±l u„=±l ^ ^ ri r„_i ™=-< 



13 

xe ^ ° Zn[C]5j2..,m., (29) 

i 

where the terms 

Z„[C]=e-^/^(^)^(^'^')^(^')'^^'^^' (30) 

are set by Gaussian fluctuations around n-instanton trajectories ©(r) ( |27|) . The 
correlator G{t,t') = {69{t)69{t')) can easily be evaluated for a dilute instanton gas 
provided both times r and r' are outside the instanton cores, i.e. |t — rj |, \t' — Tj \ ^ 
for every j. In this case Zn[C] reduces to the generating functional for a harmonic 
oscillator Zo[(] defined by Eq. ( I30l) with 



Fluctuations of PC 



8 




Figure 2. Persistent current / (measured in units eA/27r) as a function of the magnetic 
flux (j)^ for K = 3 at different temperatures: T = Q,T^ 0.125 A, T = 0.5A and T = 2A. 



where the last expression remains valid for f3Q ^ 1 and |r — r'| ^ p. Proceeding 
analogously to Ref. [H] and employing the Poisson's resummation formula we obtain 

Z[C] = Zo[C]Ee V A (32) 

fc=i 

In the limit C ~^ Eq. f l32|) reduces to the partition function [11] and the average value 
of PC is obtained from Eqs. 1^ and with = 1: 



eA ELiSin(g^)e 



/3Acos(M^i^^ 



27r(0j:-fc) - 



(33) 



At low temperatures T ^ A/k^ Eq. fl33|l reduces to a simple formula 

j=£^sin(^V -l/2<0.<l/2, (34) 

K \ K J 

i.e. at T = the magnitude of PC is proportional to A while its flux dependence 
deviates from a simple sinusoidal form for all k > 1 and tends to the saw-tooth one in 
the limit of large k. With increasing T PC amplitude decreases and the current-flux 
dependence gradually approaches the function I{(px) oc sin(27r0^). These dependencies 
are depicted in Fig. 2. 

The same equations for N = 2 yield the second current cumulant 

n(r) =V + Voscir) + Vosc{/3 - r) + P(/(r) + f{/3 - r)), (35) 

where for < t < f3 — we find 



g2A2 2. sm- ( - - J e^^'^"'^' 



A;=l 

e^A^ |'Esin(^)e^^-<^)\' 



V A:=l / 



(36) 
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and 



e^Afcti V « ) 



flAcosf^ZLWa: 



P = . (37) 

A;=l 

In Eq. (J35i) for r ^ we also defined 

/(r) = — J e{n - r)e{n)dn ^ ^e"^^ (38) 



-oo 



and 

-Poscir) = ^^^e-^\ (39) 

As one could expect from our general analysis in terms of the exact eigenstates the 
result ( 135|) indeed consists of two different - time-independent and time-dependent - 
contributions. The meaning of each of these terms can be identified with the aid of Eqs. 
(l5])-(l7l). As we already discussed, exactly at T = the time-independent part V should 
vanish, V = 0. This fact is indeed directly observed from our result (l36l) in the limit 
T 0. 

We also note that the expression for V ( l36l) can be established within the effective 
tight-binding model in which case the particle successively hops between k nods on 
the ring. This observation demonstrates that the term V is universal meaning that 
it depends only on the tunneling amplitude A but not on the profile of the periodic 
potential. Within the tight-binding model the current operator commutes with the 
total Hamiltonian, the current-current correlator does not depend on r and vanishes in 
the limit T — )■ in accordance with our general considerations. 

At non-zero temperatures, however, the term V does not vanish. AtA/K^T^f2 
from Eq. fl36l) we get 

^ = ^(/o(/3A)-/2(/3A)) 

- 2^-1 eos(2vr0.), (40) 

where Ii{x) are the Bessel functions. 

Let us now turn to time- dependent contribution to n(r). With the aid of Eqs. ([5]), 
([7]) and fl35l) we identify 

Vosc + Vfir) = ^ E |(m|/|n)|V^^"e-^(^™-^"). (41) 

The form of this - non-universal - contribution cannot be recovered within the tight- 
binding model as it explicitly depends on the instanton solution and, hence, on the 
particular shape of the periodic potential. 

In order to proceed let us observe that there exist k low-lying quantum levels 
with energies Q/2 — A cos (2-7r(03, — k)/K) in our problem. These states originate from 
tunneling depletion of the ground state energy level f2/2 in each of k potential wells. 
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Below we will label these states as |0A;) with k = 1, k. Due to the rotation symmetry of 
our model all matrix elements between these states vanish and, hence, do not contribute 
to the current-current correlation function. 

Next K energy levels with / = 1, k occur due to depletion of the first excited 
state 3fl/2. These states are characterized by the energies 3f2/2 + A cos (27r(0a; — 
where the parameter A is to be defined below. With the aid of the symmetry arguments 
one can again demonstrate that the matrix elements of the current operator between 
the states with different / vanish while the matrix elements between the states \0k) 
and remain non-zero provided k = I. In order to evaluate these matrix elements it 
suffices to consider the instanton contribution small by setting (A + A)r ^ 1 in Eq. (HTil 
and to expand the right-hand side of this equation in powers of Ar and At. Comparing 
the first two terms of this expansion with Vosc and Vf ( !37l) -( !39l) . in the limit Uq ^ Q 
considered here we identify 

Note that exactly the same expressions can also be recovered from the WKB analysis 
of the Schrodinger equation for the cosine potential. 

6. PC noise power 

Finally let us evaluate the real time current noise power defined in Eqs. (jH]), (ITU]) . 
Employing the above results at T -C O we obtain 

Atck^UqZ 

^{6{u-n-tk) + 6{u + n + tk)) (43) 
where we defined 

32^oA ( 2'n{<P^-k) \ 

^^ = ^^^°H )- ^ ^ 

We observe that - in agreement with our general analysis - PC noise power has the form 
of peaks at frequencies equal to the distance between the energy levels with non-zero 
matrix elements of the current operator plus an additional peak at zero frequency. In 
the zero temperature limit T — )■ the amplitude of this peak tends to zero along with 
the terms related to transitions to higher energy levels and Eq. fj43l) reduces to 

= ('^(^ -n-t^) + 5{uj + n + eo)) , (45) 

where eo(</'x) = niaxfcefc(0a;). This result demonstrates again that PC fluctuations indeed 
persist down to T = in which case peaks of PC noise power Si_^ occur at frequencies 
corresponding to transitions between the two lowest energy levels for which the matrix 
elements of the current operator differ from zero. We also note that S^j differs from zero 
even at zero external flux 0^. = when the average PC value is zero. In the presence of 
dissipation due to interaction of the particle with other (quantum) degrees of freedom 
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Figure 3. Zero temperature PC noise spectrum (arbitrary units) and its derivative 
with respect to the flux dSi^/d4> (arbitrary units) as functions of (measured in units 
of l/2Afi?^) for an ensemble of rings (or for a ring containing many independent 
channels) with Uq uniformly distributed within the interval from 30 /MR^ to 65 /MR^. 



the energy levels acquire a finite width, the peaks get broadened and the noise power 
should differ from zero also in a wider range of frequencies. 

Similarly, broadening of such peaks inevitably occurs in ensembles of rings or 
individual rings with many conducting channels. Within our model this broadening 
can be illustrated by considering an ensemble of rings with the parameter Uo uniformly 
distributed within some energy interval, e.g., as it is indicated in figure caption to Fig. 
3. In this case the total PC noise produced by the system is given by the sum of a large 
number of very close peaks (H5|) effectively resulting in a much smoother and broader 
noise spectrum, as it is shown in Fig. 3. 

A specific feature of PC noise is the dependence of S^^ on the external magnetic flux 
0a:. This dependence occurs due to the presence of quantum coherence in the system 
and disappears if this coherence gets destroyed. Hence, such sensitivity of PC noise 
spectrum to the flux can be used as a measure of quantum coherence in our system. 
Taking the derivative of S^j with respect to the flux, for the model considered here we 
obtain 

dS^/d(j)^ oc sin(27r0a,/K). (46) 

Typical dependence of dS^j/dcpx on u is illustrated in Fig. 3. We believe that the main 
qualitative features of our results displayed in Fig. 3 should survive also in other models 
and can be detected in experiments with nanorings. 

It is also interesting to point out a direct physical analogy between our results 
and those of Refs. fl3\ [HI |T5] where equilibrium supercurrent noise in point contacts 
between superconductors was investigated. Also in that case the noise power spectrum 
has the form of peaks which occur both at zero frequency and at frequencies equal to the 
distance between Andreev levels inside the contact. At T — )■ the zero frequency peak 
disappears while the other peaks do not vanish except in the limit of fully transparent 
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barriers. In the case of many channel diffusive contacts the supercurrent noise spectrum 
gets broadened [15] in a quahtatively similar way to the result displayed in Fig. 3. In 
addition, the noise spectrum [131 CH CH] turns out to depend on the phase difference 
across the superconducting weak link. This dependence has the same physical origin as 
the flux dependence of PC noise considered here. 

In summary, we investigated equilibrium fluctuations of persistent current in 
nanorings and demonstrated that these fluctuations do not vanish even at T = 
provided the current operator does not commute with the total Hamiltonian of the 
problem. A specific feature of PC noise is its quantum coherent nature implying that 
the noise spectrum can be tuned by an external magnetic flux inside the ring. We believe 
that the key features captured by our analysis will survive also in other models and can 
be verified in future experiments. Our further analysis will be devoted to the effect of 
dissipation on PC fluctuations in systems with many degrees of freedom. 
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